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A simple way to generate propulsion at low Reynolds number is to periodically oscillate a passive flexible 
filament. Here we present a macroscopic experimental investigation of such a propulsive mechanism. A robotic 
swimmer is constructed and both tail shape and propulsive force are measured. Filament characteristics and the 
actuation are varied and resulting data are quantitatively compared with existing linear and nonlinear theories. 
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At small scales, the physics of swimming is fundamentally different than at mesoscopic scales as the dominance of viscous 
forces over inertial forces leads to equations of motion that are time-reversible. In his famous lecture. Life at Low Reynolds Num- 
bers, Purcell Ql described three simple swimming mechanisms that are not time-reversible and hence lead to a net translation in 
the absence of inertial effects: (1) the "corkscrew" in which a rigid helical filament is rotated in a viscous liquid, analogous to 
the swimming mechanism of many bacteria ISlM; (2) the "three-link swimmer," the simplest rigid-linked mechanism that swims 
without inertia Q]; and (3) the "flexible oar" flQ-Sl, in which a flexible tail is oscillated in a viscous fluid, generating traveling 
waves along the filament that produce a propulsive force (see also The purpose of this letter is to experimentally 

investigate the flexible oar design and to compare the resulting force data with existing theories. 

Swimming at micro-scales has long been the realm of bacteria and other microorganisms but contemporary advances 

have allowed engineers to catch up with nature. Dreyfus et ai. have recently created the first manmade micro-swimmer lisj . in 
which a chain of paramagnetic beads propagates a bending wave along the chain driven by an external magnetic field. Although 
construction of this remarkable swimmer was at least partially motivated by existing flexible tail theories laRlslQ IlolfnI . the 
mechanism is not a truly passive flexible tail as internal torques are applied along the length of the filament. A second experiment 
erformed by Wiggins et aJ. measured the shape changes of a passive actin filament, oscillated at one end via optical tweezers 
. The shapes recorded in these trials match elastohydrodynamic theory well however the resulting propulsive force - a key 
parameter in designing microscopic swimmers - was not measured. Here we propose the first experimental determination of 
this force and show that the linear theory due to Wiggins and Goldstein |Q] quantitatively predicts both the shape of the elastic 
filament and the resulting propulsive, viscous forces. 

In order to experimentally quantify the propulsive characteristics of the flexible oar design, we built a robotic swimmer dubbed 
"RoboChlam" (after the algae Chlamydomonas), as is displayed in Fig.^a). The RoboChlam body was approximately 8 cm in 
length and housed a geared DC motor The motor's rotation was converted into an angular oscillation using a Scotch yoke and a 
lever (see Fig.^^). Consequently, the tail was angularly-actuated: the base of the filament was fixed at the origin and the base- 
angle was varied sinusoidally with an amplitude ao and a frequency lo. The voltage across the motor governed the oscillation 
frequency (between 5 and 0.4 rad/s), and the length of the lever controlled the amplitude of oscillation (0.814 rad and 0.435 rad). 
At the end of the lever, stainless steel wires of length 18 cm to 30 cm acted as elastic tails. Two different tail diameters were 
used in these experiments: D = 0.5 mm and 0.61 mm resulting in bending stiffnesses of 6.1 x 10^^ and 1.3 x lO^'' N • m^, 
respectively. 

RoboChlam was immersed in high viscosity (3.18 Pa • s) silicone oil to approach the low Reynolds numbers (10~^ to 10"'^) 
achieved by microorganisms. Tail shapes generated by RoboChlam were imaged with a video camera at 30 frames per second 
and 720 x 480 pixels per frame. A cantilever beam anchored the device, and a pair of strain gauges on opposite sides of the 
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FIG. 1: (a) Experimental setup to measure tail shapes and propulsive forces, (b) Scotch yoke and lever mechanism. The rotor and follower 
form the Scotch yoke, which converts the motor's rotation into a translational oscillation. This oscillation is then converted to an angular 
oscillation by a lever. The angular oscillation is approximately sinusoidal for a constant motor rotation. 



beam measured beam deflection. Strain gauge readings were converted into force measurements; a no-load voltage reading was 
taken at the beginning and end of each trial to measure the thermal drift in the strain gauges and the accompanying circuitry. 
Although force data were obtained through measured deflections of the cantilever beam, this deflection was small - less than 
half a centimeter at the beam's tip - thus, RoboChlam's position was approximately fixed. Experiments showed that an angular 
oscillation starting with the tail at rest reached steady-state motion after approximately two periods of oscillation; the time scale 
associated with this decay of transients corresponds well with the transient time scales observed in our nonlinear simulations. 
Finally, videos of the tail shapes were digitized for comparison to simulations and theoretical predictions. Experimental data are 
summarized in Figs.|2]and|3] 

To perform a quantitative comparison of experimental force and shape data with theoretical predictions, we first briefly review 
a few key results of the theory of actuated elastic filaments in Stokes flow 1?^. Consider an elastic, cylindrical rod whose 
base is attached to a fixed body (see Fig. 0J. In a low Reynolds number regime, the inertia of the fluid can be neglected and 
the fluid dynamics is well-described by Stokes equations. If the length of the tail, L, is much greater than its diameter, D, 
the hydrody namics can be further simplified by using slender body theory, the approximation of which is resistive force theory 
II12LI14II15I1 . Thus, the drag forces on the tail are linearly related to the velocity through the transverse and axial drag coefficients, 
and , respectively, and the drag force per unit length of the rod can be expressed as 

frf = -[e±nn + e||tt] (1) 

where the subscript t denotes a derivative in time, r is the position vector of a point along the tail, and n and t are the unit normal 
and tangent to the filament, respectively. We consider in this letter a planar actuation of the rod, so that n is defined without 
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FIG. 2: Force measurements for various tail lengths, L. Oscillation frequency was varied to span a range of dimensionless lengths, L where 
the dimensionless length and force are defined in equations Q and ^9} respectively. The + symbols coiTespond to D = 0.61 mm and 
ao = 0.814 rad. All other data correspond to _D = 0.5 mm and ao = 0.435 rad. There are no free parameters in the comparison between 
experiment and theory. 



ambiguities to remain in this plane. 

The elastic forces on the rod are derived from an energy functional which includes bending energy and an inextensibility 
constraint 
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where A is the bending stiffness, k is the local curvature of the tail and A is the Lagrange multiplier enforcing inextensibility. 
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Using calculus of variation we obtain the elastic force per unit length, = —6£/6r as given by | 

f, = -{Ai'sss - V'sT)n + (AV'.sV's + rs)i (3) 

where the subscript s denotes a derivative in the coordinate along the tail axis, is the local angle (see Fig.|4}, and r can also be 
interpreted as the local tension in the tail. 

Local mechanical equilibrium along the rod, leads to a pair of coupled, nonlinear partial differential equations of motion, as 
shown in 

Ipt = -J- {Alpssss - Tipss -Tslps) (4) 
+ — [A^s'^lpss + Tslps) , 
Tss - I^tV's^ = -A^^——^^s1psss-Atpss'^- (5) 

Numerical solutions to these equations were found using a Newton-Raphson iteration and are plotted along with experimental 
data in Figs. |2l and |3] 

For small deflections {i.e. assuming ^ <^ 1 such that ijj « y^), Wiggins and Goldstein 0| have shown that the motion of the 
tail can be further simplified and is described by a linear, "hyperdiffusion" equation: 

yt « --^Vxxxx (6) 
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FIG. 3: Comparison between experiment, linear and nonlinear theories of tail shapes. Snapshots are shown at four points in the cycle for 
one tail with L = 20 cm, D — 0.5 mm, ao — 0.435 rad, at three different oscillation frequencies: (a) ui = 0.50rad/s (£ = 1.73), (b) 
uj = 1.31 rad/s (£ = 2.20), (c) uj = 5.24rad/s (£ = 3.11). 



where subscripts x and t denote derivatives in position and time, respectively. For the case of harmonic angular-actuation, we 
apply the boundary condition ip = uq sm{ujt) at the base. The nondimensionaUzation of Eq. (|6j is obtained by substituting 
X = Lx, y = aoLy, and t = i/u into Eq. leading to yi « ~{£uj/ L)^yxx£s where, £^ = {A/uj^^y^^, is the characteristic 
penetration length of the elastohydrodynamic problem; solutions to Eq. decay exponentially in space over this typical length 
scale. The time-evolution of the tail shapes is then only a function of the angular amplitude, ao, and the dimensionless length, 

t \ 1/4 



C = L/C = L i^-^j . (7) 

This dimensionless length is the key parameter in the problem and represents the "floppiness" of the tail and hence the overall 
effectiveness of the swimmer In particular, theory predicts an optimal dimensionless tail length as both short, stiff tails and long, 
flexible tails produce negligible net translation - the first is ineffective owing to the scallop theorem and the second owing to the 
excessive drag on the long passive filament. 

For a tail that is periodically oscillated, Eq. (|6} can be solved analytically At the base of the filament, the reaction forces 

and torque must balance the drag forces along the tail. The opposite end of the tail is force and torque-free such that tps — 0, 
^ss = 0, and r = at ,s = i. For small deflections, the x-component of local drag force, Eq. Q, can be integrated along the 
length of the tail to yield the propulsive force 

{F) « -A^^—^{y,y,,,, - iyxx')x=0 (8) 
4-L ^ 



where (...) denotes averaging over one period of oscillation. Note that Eq. differs from the one presented in llaLZl] by a factor 
(^j^ — ^11 this disparity arises from a proper integration of the drag force on the filament 0]. 

In order to obtain results valid beyond the small-slope approximation, our numerical solutions to the full nonlinear system 
(Eqs. 0]and|5} were employed and the propulsive force was found by numerical integration of the local drag force along the 
length of the tail. A dimensionless force T was defined by substituting x — l^x, y = aoCy, and t — 2ni/u} into Eq. (|8}, such 
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that 

{F)=a,HJ{U-(^^)\iu\{T). (9) 

Since the distance from the tail to the nearest wall was on the order of the tail's length, drag coefficients corrected for wall effects 
as in llisll were used in simulations and for nondimensionaUzing the force data. These wall corrections have been shown to 
match well with experimental results ; for simplicity, the effect of only a single side-wall was considered. These equations 
produced a drag difference of approximately ^± — ^|| = 3.35 Pa ■ s - about 40% greater than the drag difference without wall 
effects. 

The results of our investigations are summarized in Figs.|2Eland|5] We first display in Fig.|5]the propulsive force generated 
for a range of dimensionless tail lengths, C All parameters of the experiment were known or measured, and no fitting of data 
was necessary. We obtain excellent agreement of the propulsive force with the theoretical (linear model, Eq. |6j and numerical 
values (nonlinear model, Eqs.|3and|5}- The force data from the RoboChlam experiments show a maximum dimensionless force 
at £ « 2.1, in agreement with prediction from the theory. Note that our data was nondimensionalized with the drag difference, 
£,± ~ (see Eq.|9}, instead of the transverse drag which was used in The drag difference orginated in Eq. (|8}, and it 

represents the correct scaling as a tail with isotropic drag (^j^ ^ ^||) should produce zero propulsive force 01^- We note also 
that the maximum value of C that could be tested was limited by the motor's rotation rate and the length of tail that would fit in 
the experimental apparatus. 

In comparing the data to linear elastohydrodynamic theories, there are three primary sources of error: wall effects, thermal 
drift in the experiment, and the neglected nonlinearities in the theory. The error bars in Fig.|2]arise from uncertainty in the no- 
load voltage of the strain gauge measurements. At lower oscillation frequencies, the sample time of the experiment increased, 
leading to noticeable thermal drift in strain gauge (force) measurements and thus, larger drift error for the left-most points of a 
given data-set. In addition, the tip of the longest tail (30 cm. A) was only a few centimeters from the back wall and thus, this 
wall had a non-negligible effect on the drag of the longest tail resulting in an increased thrust as expected. Recall that our wall- 
corrected drag coefficients only account for a single wall - the side wall rather than the back wall - of the tank, as is appropriate 
for all but the longest tails in our experiments. It is interesting to note that, in these experiments, nonlinear effects are completely 
negligible relative to the other two sources of error even for long tails and large actuation angles. 

In Fig.l^we plot the tail shapes from experiments along with simulationed shapes from both the linear and nonlinear theories. 
The plot shows three tails from a single data set (constant L, D, and oq, but varying uj) with dimensionless lengths (a) C = 1.73, 
(b) £ = 2.20, and (c) £ = 3.11. These dimensionless lengths span the region near the maximum dimensionless force. The 
tail shapes from experiment matched well with those from the linear and nonlinear simulations, and only slight differences 
between the three tails were observed. Tails whose dimensionless length was small (Fig. |3^) moved stiffly, while those with 
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FIG. 5: Normalized, time-averaged differences between linear (L), nonlinear (N), and experimental (E) tail shapes. The difference is calculated 
as the ^^-norm of the vertical-distance vector between two tails divided by the tail length and the number of points along the tail. Two data 
sets are shown: (1) i = 20 cm, D = 0.5 mm, and oq — 0.435 rad; (2) L = 18 cm, D = 0.63 mm, and — 0.814 rad. 



large dimensionless lengths (Fig. |3j;) were flexible, as predicted by theory. The difference between the different tail shapes 
(theory, experiments, simulations) is quantified in Fig. |5] The measured errors are observed to be small. The fact that the data 
match the linear simulation better than the nonlinear solution is fortuitous and merely reflects the fact that resistive force theory 

n 

is only an approximation of the equation of hydrodynamics 112]. 

In summary, we have presented an experimental investigation of Purcell's flexible oar swimmer Measurements of propulsive 
forces and time-varying shapes are in agreement with the results of resistive-force theory. Remarkably, the small-slope model of 
Wiggins and Goldstein Q] appears to remain quantitatively correct well beyond its regime of strict validity. 

Our future work will investigate the efficiency of this propulsive mechanism when embedded in a synthetic free-swimmmer - 
that is, an elastic filament attached to a body which translates and rotates with the forces and torque generated by the propulsive 
tail. Preliminary free swimming experiments show that rotation of the swimmer body significantly changes the shapes of the 
tail, modifying the force curve shown in Fig.|2j and appreciably impacting the dynamics of the swimmer. 
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